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Abstract 

In this paper a new method for computation of higher order cor- 
rections to the saddle point approximation of the Feynman path in- 
tegral is introduced. The saddle point approximation leads to local 
Schrodinger problems around classical orbits. Especially, the saddle 
point approximation leads to Schrodinger problems around classical 
periodic orbits when it is applied to the trace of Green's function. 
These local Schrodinger problems, in semiclassical approximation, can 
be solved exactly on the basis of local analytic functions. Then the 
corrections of the semiclassical result can be treated perturbatively. 
The strength of the perturbation is proportional to h. The pertur- 
bation problem leads to ordinary differential equations. We propose 
these equations for numerical calculation of corrections, since they can 
easily be solved by computers. We give quantum mechanical general- 
izations of the semiclassical zeta functions, spectral determinant and 
trace formula. 



Feynman's path integral in quantum mechanics and similar path inte- 
grals for stochastic systems are the most intuitive tools of modern theoretical 
physics. But calculations with path integral are difficultly in general. One 
can often find numerical solution of the underlying partial differential equa- 
tion easier. The most convenient asymptotic method to evaluate the path 
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integral is the saddle point approximation. The leading gaussian approxima- 
tion is easy to perform and gives very good results. This is the case when 
we calculate the trace of the Green's function of the one-body Schrodinger 
equation in gaussian approximation and get the Gutzwiller trace formula [Q]. 
From the Gutzwiller trace formula one can recover the eigenenergies of the 
Schrodinger equation for bound systems and the resonances of the S ma- 
trix for scattering systems in semiclassical approximation]^. There are many 
attempts to improve the semiclassical approximation to get more accurate 
energies and resonances. But the gaussian approximation has its inherent 
limitation and one should go beyond it to improve the accuracy. Recently, 
Gaspard et al.]^ computed corrections of the Gutzwiller trace formula and 
showed that the resonances of the two and three disk scattering systems |TT| 
can be improved considerably. They have used the usual Feynman graph 
technique of the perturbation theory and computed large number of graphs 
to get the corrections. In general the conventional graph calculus is very cum- 
bersome. We give here two simple examples, where the corrections can be 
computed and an agreement with the calculations in Ref. can be checked. 

In this paper we show how the corrections to the gaussian approximation 
can be computed using simple ordinary differential equations. We hope the 
results can be generalized for the field theory and for many body systems 
in the near future. Here we treat the one body non-relativistic Schrodinger 
equation and give the h corrections of the Gutzwiller trace formula and the 
spectral determinant introduced by Voros|0]. 

1 Saddle point evaluation of the trace 

The time domain propagator G{q, q', t) of the Schrodinger equation is defined 
by 

ihdtG{q, q', t) - HG{q, q', t) = 5{q - q')5{t), (1) 

where 

H = -^A + U{q) 

is the Hamilton operator in units of m = 1. The path integral representation 
of the propagator 

G'(g,g',t) = y'Pg"e^^('?''?'*'^"), (2) 
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where / T>q" represents the functional integral measure for all the paths con- 
necting q with q' in time t, and S{q, g', t\q") is the classical action between q 
and q' computed along a given path q" . The resolvent of the Hamilton op- 
erator, the energy domain Green's function, can be recovered as the Fourier 
transform of the propagator 

1 fx \ 

G{q,q',E) = -j^ dtexpi^-EtjG{q,q',t). (3) 

The trace of time and energy domain Green's functions can be expressed 
with the help of the eigenenergies of the Schrodinger equation: 

TrG(t) = f;exp(-^E„t), (4) 



n=0 

OD -j^ 



n=0 E-^n 

The path integral expression for the trace of the propagator is 

TrG{t) = J dqG{q, q,t) ^ J Pg"et^(*l''"), (6) 

where / Vq" now represents the functional integration for closed paths. The 
spectral determinant A(£') = det(£'— ^) = UniE—En), whose zeroes are the 
eigenenergies or resonances, is related to the trace of energy domain Green's 
function through the logarithmic derivative 

TTG{E) = ^logA{E). (7) 

In the saddle point approximation we assume that the leading contribution 
to the path integral is coming from the neighbourhood of paths for which 
the classical action is stationary. This condition singles out the classical 
trajectories from the infinite variety of possible paths. The trace of the 
propagator can be written as a sum for the classical periodic orbits: 

TrG(i) = E / ^?pexp (^^W,{x,t)j , (8) 

where Z^p.o. denotes the summation for the classical primitive periodic orbits, 
/ T>qp denotes a functional integral in the neighborhood of periodic orbits 
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and Wp{x,t) = Y.n'u^n{x — q{t)Y /n\ denotes the classical action expanded in 
multi-dimensional power series around the periodic orbit. The energy domain 
Green's function can be recovered by Fourier transformation. 

Here we do not investigate the general validity of saddle point approxi- 
mation. However, it is important to note that the regions around periodic 
orbits, where the action is expanded in power series, may overlap and this 
causes some overcounting in the formula (§). Therefore in computations the 
number of periodic orbits included in the sum should depend on the order of 
truncation of the power series. In the semiclassical or gaussian approximation 
the method proposed by Berry and Keating[^ can deal with this problem. 
We hope, that their method can be extended for the situation discussed here. 

Now, the traditional computation goes as follows: In the gaussian ap- 
proximation the power series expansion is truncated at quadratic terms and 
the Gauss integrals are evaluated. This way one can derive Gutzwiller's trace 
formula. Corrections of the gaussian approximation can be found by expand- 
ing the action to higher orders, expanding the exponential and performing 
the gaussian cumulant integrals. 

However, we have a better choice at this point. We can realize that the 
terms in (^) are similar to the original one in the neighborhood of periodic 
orbits except the classical actions are given in power series form. We can 
look for the partial differential equations corresponding to these local path 
integrals. These partial differential equations are Schrodinger equations in 
a power series form around periodic orbits in the configuration space. The 
terms in (||) are traces of the propagators TrGp(t) of these local Schrodinger 
problems 

TrG'(t) = TrGp(t). (9) 

P.O. 

By Fourier transformation we can also get the energy domain trace term by 
term 

i:iG{E) = J2 TrGp(£;). (10) 
p.o. 

The Gutzwiller-Voros spectral determinant is then a product of the spectral 
determinants of local Schodinger problems 

A{E) = n \iE), (11) 
p.o. 
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which are related to the traces as logarithmic derivatives 



TrG,(E) = ^logA,(E). (12) 

The product (|TT|) is the direct generalization of the Selberg type product [|, [7| 
as we will see later. 

In general, there are many different types of closed periodic orbits which 
can contribute to the product (]TT|). Some of them are zero length or equ- 
ilibrium[Q orbits which can be well treated by conventional methods. The 
spectral determinant of the zero length orbits gives a smooth contribution, 
which is the counterpart of the Weyl or Thomas-Fermi terms. From now on 
we neglect these terms since they do not change the location of the zeroes 
of the spectral determinant. Also, the periodic orbits in the complexified 
phase space of the hamiltonian system can occur Special contributions 



can come from diffraction cycles introduced recently [13, 11, |T^, 18 1. Here 
we concentrate on usual classical periodic orbits and we discuss other terms 
elsewhere. 



2 Local spectra of the Schrodinger equation 

To compute the local spectral determinants Ap{E) we have to solve the local 
Schrodinger problem in the neighborhood of a classical periodic orbit. The 
most convenient way to do this is to rewrite the Schrodinger equation 

ihdtip = Alp + Uip (13) 

2m 

with the usual ansatz 

^ = cpe'^/'', (14) 
which yields the following equation 

- fdtS + ihdtip = (a^ + 2i/hVipVS + i/tupAS - i/n^^{vsf) + U^. 

(15) 

Here we have many possibilities to group the terms since we have not made 
any restriction for S and (p. Our main concern is to separate the classical 
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and the quantum time evolution. Therefore, we require the phase to fulfill 
the Hamilton- Jacobi equation 

dtS + ^{VSy + U = 0. (16) 

At a given phase S{x, t) the complex amplitude fulfills 

dtip + VipVS + ^ipAS - yAv? = 0. (17) 

This is a Fokker-Planck like equation in the velocity field of the classical 
action V5'(x,t) with a complex diffusion constant ihj2. 

If the local Schrodinger equation around the periodic orbit has an eigenen- 
ergy E the corresponding eigenfunction fulfills 

^p(t + Tp) = e-'^^-/Vp(i)- (18) 
For a general energy value one can find eigenf unctions '?/'p(t) such that 

^J(t + T,) = e-^^V^Ai(£;)V'J(t). (19) 

If the eigenvalues Ap(i?) are known the local functional determinant can be 
written as 

A,{E) = Yi{l - \\{E)\ (20) 

since Ap(i?) = yields the eigenenergies of the local Schrodinger problem. 
We can insert the ansatz (|1^) and reformulate (^) as 

em^^v)^i^[^t + Tp) = e-^'^^-l''\\[E)e^^^^^\{t). (21) 

The phase change is given by the action integral for one period S{t + Tp) — 
S{t) = Jo^ L(t)dt. Using this and the identity for the reduced action Sp{E) 
of the periodic orbit 

Sp{E) = jpdq = J^" L{t)dt + ETp, (22) 

we get 



e^^^(^VUt + T,) = A;,(E)^;,(t). (23) 
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Introducing the eigenequation for the amphtude 

cplit + T,)=RliE)cp^^it), (24) 
the local spectral determinant can be expressed as 

A,{E)=llil-RliE)e^nSAE)y (25) 

At this point we can reexpress the Quantum Gutzwiller-Voros spectral 
determinant with the local eigenvalues, which now reads as 

m) = I[W-K(E)e'''''^^^). (26) 

p I 

The trace formula can be recovered from (|^: 

1 i?i(E)e^^''(^) / , , d\ogRi(E)\ , ^ 

We can also rewrite the denominator as a sum of a geometric series and get 
TrG(i?) = ^EE (^^(^) -^^^^^^^1^) E(^p(^))^e^^^(^). (28) 

The new index r can be interpreted as the repetition number of the primitive 
orbit. This expression is the generalization of the Gutzwiller trace formula 
for the exact quantum mechanics. We can see, that the h corrections enter 
here in a more complicated form than in the spectral determinant. 

The local eigenvalue spectra can be computed with a variety of numer- 
ical methods. These numerical eigenvalues then can be used to build up 
the spectral determinant or the generalized trace formula. For theoretical 
calculations we develop here an analytic perturbation method. 



3 Expansion in h 

Since ^ is a small parameter we can develop a perturbation series for the 
amplitudes 

^'W = E T v'^'^Ht) (29) 
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which can be inserted into the equation (|T7|). The eigenvalue can also be 
expanded in powers of ih/2: 

Ri{E) = e^v[t^{^fj^ Ct^y (30) 
The h expansion of the eigenvalues is 



m 



ME) = explf (|) Cr^\ (31) 



.m=0 



= exp(Cr)(l + |c,'"+(f)'(i(C;")^ + cP')+.... (32) 

The eigenvalue equation in h expanded form reads as 

£ (|)%'M(, + T,) = exp I £ (jYct'] . E {j)\''-\t). 

(33) 

Expanding the eigenvalue like in (|32D and collecting the terms of the same 
order in h yield a set of eigenequations 

^'W(t + T,) = exp(Cj°V'^°Hi), 

^'«(t + r,) = exp(c;°))b'«(t) + c(\'(°)(t)], 

+ T,) = eMCPW^'\t) + C/\'«(t) + (Cf) + i(Cp))^)^'W((fi,) 

and so on. These equations are the conditions selecting the eigenvectors and 
eigenvalues and they hold for all t. Without loss of generality we can also 
assume that ip''^'^\0) = 1 and (^"^'"-'(O) = for m > 0. By adding these 
assumptions we can simplify the equations 



^'W(T,) = expiCP), (35) 
^'(^)(T,) = eMCl'^)Cl'\ (36) 
^'(^)(T,) = eMCP)iCl'' + liCiy). (37) 
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4 The analytic eigenbasis 



In the neighborhood of a classical periodic orbit we can look for the solution 
of the Hamilton- Jacobi equation in a power series form. Let q{t) denote a 
classical periodic orbit with period Tp. Let us expand the phase around the 
time dependent trajectory: 

oo 

Six,t) = J2s^{t){x-qit)r/n\, (38) 

n 

where s„ are d"' matrices in dimension d with time dependent entries. This 
power series ansatz assures that the local Schrodinger problem corresponds 
to the one in (|). To derive ordinary differential equations for the expansion 
coefficients let us expand the potential around the periodic orbit 

oo 

Uix)=Y.Ur,{t){x-q{t)r/n\, (39) 

n 

where d^ matrices in general. If we put this two expressions in the 

Hamilton- Jacobi equation we get in the one dimensional case 

1 " n\ 

Sn - Sn+iq + T^Yl r~rrn7l +Un = 0. (40) 

In the multidimensional case we get similar expressions for the entries of the 
s matrices. If the si vector is choosen to be the momentum of the classical 
orbit 

p = q = si, (41) 

the equations are simpler and their meanings are obvious. The first eq. in 
the hierarchy corresponds to the classical action along the path: 

2 

So = j-Uo = L{t), (42) 

where L{t) is the Lagrange function evaluated on the periodic orbit. The 
second is the Newton equation 

P=-ui, (43) 
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since -ui is the force along the trajectory. The dxd matrix S2 is famihar from 
the wave packet theory and describes the shape of a gaussian wave packet [p^] 



S2 = -sl - U2. (44) 

Trs2 describes the expansion of infinitesimal volume elements evolving along 
the classical orbit. The next equation is 

S3 = -3S2S3 - Ms (45) 

and the rest of the equations are linear equations for s„. These are pure 
classical equations describing the analytic structure of the action around 
periodic orbits. 

The gradient of S{x,t) entering in the amplitude equation (|T7p are pe- 
riodic along the periodic orbit. Therefore, the s„, > matrices are also 
periodic yielding the boundary conditions Sn(t) = Sn(t + Tp) where Tp is 
the period of the orbit. The term so(t) is not periodic and is given by the 
action integral so{t) = Jq L{t)dt. si{t) is the momentum along the periodic 
orbit and it is a periodic function. The most complicated equation is (|4^). 
In general it has more than one periodic solutions. In case of unstable pe- 
riodic orbits the solution of the equation (^) converges to a single stable 
solution starting from almost all initial conditions. The rest of the solutions 
are unstable. The wave packet described by the stable solution is decaying 
in time, while the rest of the solutions describe wave packets with increasing 
amplitudes. These solutions are non-physical, since they describe local wave 
functions with exponentially increasing norms. We have to exclude these 
solutions. In case of stable periodic orbits we also have only one solution 



of (S) for which the local wave function is decaying and we have to choose 
this solution. The higher order {n > 2) equations are linear in s„ and their 
unique periodic solutions can be found order by order. 



5 Evolution of the amplitude 

After solving locally the Hamilton- Jacobi equation we can look for the local 
solution of the amplitude equation. We can expand the amplitude around the 
classical path in power series. This analytic basis is appropriate for classical 
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Perron-Frobenius operators since it is very easy to diagonalize the evolution 
operator on this basis Inserting the expansion 

oo 

^{x,t) = J2an{t){x-q{t)r/n\ (46) 

n 

into the equation (|T^) yields the following equations for the coefficients in 
one dimension 

n\ / 1 \ ih 

an - an+iq + Y] 7 77777 Ctn-i+lSz+l + -an-iSi+2 - — = 0. (47) 

J^{n- ly.ll V 2/2 

In the multidimensional case we get similar equations for the expansion co- 
efficient matrices. Using eq. (|4l|) one can shghtly reduce these equations: 



^2 , iJi 

ao = - Y«o + ya2, (48) 

3s2 S3 ih 
oi = — ^ai - — ao + —as, (49) 

5s2 S4 ih , . 

02 = — —a2 - Zssai - —ao + —a^, (50) 

and so on. These equations are linear and have the general form 

(2n+l)s2 ih 
^ -an--- + iran+2- (51) 



2 2 

These equations are matrix equations in higher dimensions but their structure 
remain unchanged. 



6 Stationary solutions 

The set of equations (^D]),(^^ and (|T^) define the full set of equations we have 
to solve. This can be carried out without further considerations. However, we 
know beforehand, that the solutions of the Scrodinger equation are stationary 
states. The stationarity condition implies that the phase of the wave function 
fulfills the condition 

= -E (52) 
dt ^ ^ 



11 



and the amplitude has no exphcit time dependence 



t) 

dt 



0. 



(53) 



These equations give us some additional equations for the expansion coeffi- 
cients, which have the form 



in the one dimensional case. In the multidimensional case the coefficient 
matrices fulffil similar equations. These equations can help us to reduce the 
equations we really have to solve, since some of the higher order expansion 
coefficients can be expressed by the time derivatives of the lower order coeffi- 
cients. In one dimension all the higher coefficients can be directly computed 
from the time derivatives of the zero order terms. In two dimensions the 
number of s„ and a„ matrix elements is n + 1. The number of the additional 
equations derived above is n. Therefore, on each level we get 1 entirely new 
equation, which we have to solve. In three dimensions we get n entirely 
new equations for the phase and the amplitude on each level. We show on 
examples how the reduction of the equations can be carried out. 

7 h expansion on the analytic base 

Inserting the expansion (|29|) in the equations (^^^-|5T|) we get coupled equa- 
tions. The zeroth order or semiclassical equations form an autonomous sys- 



Sn{t) - q{t)Sn+l{t) 

anit) - q{t)an+i(t) 



for n > 0, 
for n > 0, 



(54) 
(55) 
(56) 



tem 



«n - a^+iq + 2^ _ ^yj^^ [a^-i+i'Si+i + -a^n-i^i 



:) = 0. (57) 



For example, the ffist three equations have the form 




(58) 




(59) 
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= -^af -2.3aS°)-|ar. (60) 

The important feature of these equations is that they are hnear and have the 
general form 



(0) _ (2n + l)g2 ^(o) 



2 

and so on. 



af + -, (61) 



We can solve equation (^) and get 

aT\t) = aT\0) exp (- ^s^m^ (62) 



from which one can read off the eigenvalue 



Cl,'^ = - ls,{t)dt. (63) 
Jo 2 

The S2 in general goes through 1/t type singularities. If this happens the 
integral should be carried out by principal value integration. This means 
that we add a regularizer ie term. The integral then can be calculated and 
in the e — limit we recover an i7r/2 Maslov phase. 

The rest of the equations do not play a role in yielding the first eigen- 
value. The solution aQ^^\t) can be inserted into the next equation (^9]). Since 
equation ( |5^ is a linear one driven by aQ^^\t) its particular solution fulfills 
the condition (|34D. The rest of the equations can be solved the same way and 
we get the eigenamplitudes al^'^\ The rest of the semiclassical eigenvalues 
can be recovered by setting a^^^-* = for n < I. Then the /-th eigenvalue is 
given by 

cr = f' s.m- (64) 

The semiclassical eigenvalues are connected with the stability properties 
of the periodic orbits. For example the first (/ = 0) eigenvalue is related to 
the product of the expanding eigenvalues |]T5| 



^MCS"') = JYiJ^2 (65) 

where Aj denotes the expanding (Aj > 1) eigenvalues of the linear stability 
or Jacobi matrix of the periodic orbit and Up is the Maslov index of the 
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periodic orbit. The Maslov phase comes from the singularities of S2{t) (See 
e.g. Ref[[l5l.) The product 

HE) = nn(l - exp{zS,/h + Cj'^)) (66) 
p I 

in this approximation is known as the Selberg product p, |^. 

After the calculation of local semiclassical eigenvalues and eigenvectors we 
can use them in the next level of approximation. The differential equations 
connecting the m + 1-th order amplitudes with the m-th order amplitudes 
are 

(m+1) J"^+^)a I ^' I ^-n'-"^^ 



Again, using Eq. (^) one can slightly reduce these equations 

. (m+1) _ S2 (m+1) (m) /^^jx 

— "y'^O + '^2 J iDoj 

. (m+1) _ 3S2 (m+1) S3 (m) . s 

Ol — ^fli —Uq + , {Oy) 

(m+1) _ 5S2 (m+1) 9^ ^(m+1) 54 (m+1) (m) ,„f.s 

0-2 — ^'^2 ~ ^S^tti — —% +04 , (_/UJ 



and so on. These equations are linear and have the general form 

a^^) = - (^I^±ll£^ar^) + ... + atl (71) 



Inserting the eigenamplitudes a5|™''(t) we get linear driven equations for the 
next order of the amplitudes. The solutions of these equations, which satisfy 
the conditions of type (|34D , yield the corrections Ci"^^ of the semiclassical 
eigenvalues C;''°\ 

It is important to note that the driving terms in these equations cause 
resonances, since they are solutions of the homogenious equations they drive. 
The corrections ajf™-* (t) are not exponentially growing functions of time. In- 
deed, a solution of ( |68D behaves like 

af\r ■ Tp) ~ exp(rC;°^)(a + 6 ■ r), (72) 



14 



due to the resonance, where a and b are appropriate constants. The particular 
solution satisfying (|35|) is the one which does not contain pure exponential 
part (a = 0) and then the correction can be read off as C;*-^^ = b. 

As a consequence of this hierarchy, it is increasingly difficult to get cor- 
rections of the eigenvalues corresponding to large /. It is more convenient 
to reorganize the quantum Selberg product as a product of quantum inverse 
zeta functions 

AiE)=l[Cr\E) (73) 
where the quantum zeta functions are defined by 

Cr\E) = n f 1 - exp{^S,{E)/n + J:{^m^C!^"'\E))) ■ (74) 

p \ m / 

These zeta functions are the quantum generalizations of the Ruelle zeta 



functions[114|, iLuj . The leading resonances and the eigenenergies can be com- 
puted from the zeroes of the / = quantum zeta function. The curvature 
expansion proposed by Cvitanovic and EckhardtQ can also be applied using 
the new quantum mechanical weights 

t, = exp (^S,{E)/h + Y.{^h/2rC',^"'\E)] . (75) 



8 Billiards 

We have to discuss some special features of the billiard systems here shortly. 
In billiards the potential is not an analytic function, therefore the theory 
can not be used without further considerations. We have Dirichlet boundary 
condition for the wave function on hard walls. The wave function should 
vanish on the wall. Our approach here is basically the tracing of a wave 
packet along the classical trajectory in the configuration space. When the 
packet is hitting the wall, the incoming wave function at time t is given by 
the packet right before the collision, evaluated on the wall 

i:Ux{s),y{s),t) = ^(x(.), !/(.), ^_o)e^^«^(^)'^(^)'*-»)/^ (76) 

where {x{s),y{s)) is some analytic parametrization of the wall around the 
classical point of reflection. The outgoing wave function is the wave function 
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right after the colhsion 

^^^ouMs):y{s):t) = (^(x(s), t+o)e^^«^^^)'^(^)'*+°)/". (77) 

The sum of the incoming and the outgoing wave functions should vanish on 
the hard wall due to the Dirichlet condition. This implies that the incoming 
and the outgoing amplitude and the phase has the relation 

ip{x{s),y{s),t-o) = (p{x{s),y{s),t+o) (78) 

and 

S{xis),yis),t_o) = Six{s),y{s),t+o) + ^vr. (79) 

These relations mean, that the power series with respect to s of these func- 
tions are equal on both sides of the collision modulo the tt phase shift. This 
phase shift can be interpreted as the Maslov phase coming from the hard 
wall. 



9 Example I. 

As an example first we study a simple two dimensional potential system 
with potential U{y) = —^X^y^ — -^ay*i which is bounded by hard walls in 
X — Q and X — L (Fig. 1). The only classical periodic orbit in the system is 
lying on the y = line and the particle is moving with constant momentum 
p = \/2E, (m = 1) along it. 

The equations for the y derivatives of the phase are 

Syy + Syy + Uyy — 0, 

^yyy "I" ^^yy^yyy "I" ^yyy ~ 0; 
'-'yyyy ~ ^'-'yyyy'-'yy ~ ^'~^yyy ~ yyx ~ '-'yyyy ~ 

and so on. We will compute here only the first correction, so we have to solve 
only these equations. The stationarity conditions yield 



^yy 


pSyyx 0, 


(83) 


Sy 


~ pSyx = 0, 


(84) 




pSxx — 0, 


(85) 


Sxyy 


pSxxyy 0. 


(86) 




16 





(80) 
(81) 
(82) 



The classical momentum is the first derivative of the phase 

Sx = p, (87) 
Sy = 0. (88) 

Since these components are constant in time Eq. (|8^) and ( pS]) yield S^y = 
and Sxx = 0. The periodic solutions of (|80D are constant in time and given 
by Syy = ±A. The stable solution is the one with positive sign. The periodic 
solution of (|8l|) is Syyy = 0. Since Syy is constant (|83|) gives us S^yy = 0. The 
solution of (0) then simply Syyyy = a/4A. 

The zero order amplitude equations for the y derivatives are now 

a + ^Syya = 0, (89) 
3 1 

% + 2^yy^y + 2^yyy^ = 0' (90) 

5 1 

^yy ~^ "^^yy^yy ~^ ^yyx^x ~t~ ^^yyy^y ~^ ^(*^yy^^ ~^ ^yyyy^^ 0' (^1) 

(92) 

The solution of the first one is a = exp(— At/2). In one period (Tp = T = 
2L/p) this solution decays as ci(T) = exp(— ATp/2) so Cq*^ = —XT. Similarly, 
C/°^ = — (/ + 1/2) ATp. The semiclassical spectral determinant for this system 
then reads as 

A{E) = n (1 - exp{i2pL/h - (/ + l/2)Ar)) . (93) 

The first semiclassical zeta function is the term with index / = 0: 

C,\E) = (1 - exp{z2pL/h - l/2Ar)) . (94) 

Now we compute its first quantum correction. To do this, we have to solve 
the first correction equation 

a(^^ + k.a^^^ = <^ + ai°i- (95) 



2-yy" yy 



The function ayy is the solution of (0) which decays as ayy{nT) = exp(— AnTp/2) 
in time. This solution is 



<^yy{t) = -e"'*^'^- (96) 
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The be calculated using the stationarity relations 



a — pax = 0, (97) 
dx - paxx = 0. (98) 

(99) 



The result is 

A2 



axx{t) = -e"^*^'^- (100) 



The resonant solution of the correction equation (p5| ) is then 

a«(t)=te-^*/2f^-^V (101) 

li6A vy ^ ^ 



From here we can read off the first correction as 



a A2 



The corrected first zeta function then reads 

Co\E) =(^1~ eM^2pL/h - 1/2XT + f (j^ " ^)) • (103) 

This procedure can be continued step by step. From the zeroes of the first 
zeta function (^^{Ei) = we can recover the leading resonances of the sys- 
tem, which are the closest to the real axis and have the longest lifetime. The 
appearance of a in the expression shows, that the first correction is sensitive 
for the neighborhood of the periodic orbit. 



10 Example II. 

Our second example is an open billiard system depicted on Figure 2. The 
curved wall is locally given by the symmetric curve: 

x{y) = lc,y' + ^C,y' + .... (104) 
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The separation between the two opposite tips is L. The only periodic orbit is 
the particle bouncing force and back between the tips with constant momen- 
tum. This situation occurs for example in the two disc scattering system 
and in the confocal hyperbolae problem||22|, where the first correction to the 



Gutzwiller trace formula have been computed numerically. Since the system 
is symmetric, it is more convinient to work on the left part of Fig. 2 and 
assume, that a vertical straight hard wall is placed on the symmetry axis. 
The expansion coefficients of the phase and the amplitude fulfill Eqs.(|80|- 



9T|) without potential term {U = 0). In addition to these equations we have 
to use the bounce conditions ([78| ) and (|79|). If the classical trajectory is 
starting from the curved wall, the bounce conditions give us the following set 
of equations: 

Syy{T) + C2S^,{T) = Syy{0)+C2SM, (lOS) 
SyyyiT) = Syyy{0), (l06) 
SyyyyiT) + QC^Syy^iT) + C ,{T) = S yyyy{Q) + QC yy^ + C 

<T) = ^a(O), (108) 

ayy{T)+C2a,{T) = ^{ayy{0) + C^a^) , (109) 

(110) 

where 

|A|V2 = e^o°\ 

The solutions of the amplitude and phase equations are 

Syy{t) = (111) 

Syyyit) = 0, (ll2) 

B 3 
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452(t + to)5/2' (117) 
ay{t) = 0, (118) 

""yy^^) = (t + to)5/2 + (t + to)7/2- (119) 

(120) 

The constants B, C, to, A can be set by the bounce conditions ( |105| - |109|) and 

we get 



A = C2L + 1 + ^(^2^ + 1)2-1, (121) 
_ 3L 1 2C4L4 

p3A2(l_l/A4) +p3(A- 1)4(1 _1/A4)' ( ) 

^ _pSAi-^ 1 

2L 1 + 1/A V 

(125) 

The solution of the first correction equation (P3| ) now reads 

aW(t) = [t + to)-'/' f dt\aUt') + ayy{t')){t' + tof'\ (126) 



From the solution we can read off the first correction, which is finally 

This result can be compared with the findings of ref.[^ for the two disk 
system and with those of ref. for the confocal hyperbolae. In case of the 
two disk scattering system C2 = l/a and C4 = 3/a'^, where a is the radius of 
the disk. In this case we get 

Ci'^ = (128) 
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which coincides withe the result of Ref . Q derived via Feynman graph tech- 
nique. In case of the confocal hyperbolae C4 = — 6C|/L and the correction 
is 

= ^. (129) 
P 



This was numerically confirmed in ref. P2| . 

An important consequence of the result ( |127] ) is that in the L — > limit, 
when the walls are touching each other, it gives the universal limit 

Cj^) = 5/8a. 

In this limit the particles scattering from outside spend a long time in the 
horn formed by the touching walls. This can be considered as the prototype 
of the intermittent motion. In this limit the first zeta function in leading 
order reads as 

C„-(B) = l-e.p(>A). 

The zeroes of this zeta function are 

5 



The energy spectra is given by 

1 



27r n + 1 



2 \ 167ra(n + 1 



We can expect, that in case of intermittency due to touching billiard walls, 
the local energy spectra has an universal Rydberg type character. 



11 Conclusions 

In this paper we have proposed a new method to evaluate corrections to 
the leading saddle point approximation of the Feynman path integral. The 
method reduces the problem to a set of ordinary differential equations which 
have to be solved at certain boundary conditions. In all orders the prod- 
uct structure of the functional determinant A(i?) is maintained. One can 
introduce the quantum zeta functions. The corrections to the leading zeta 
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function is easier to calculate than a general / > term and it is very practi- 
cal to use it for extended computations. A nice feature of the method is that 
the problem of marginally stable periodic orbits can also be treated, since 
the differential equations give nonzero contributions (C; ^ 0) for / > in 
case of such an orbit. The method is applicable also for general second order 
partial differential equations, where the diffusion term is weak. Since the 
individual saddle point terms TrGp(£') are more accurately calculated than 
in the Gauss approximation one can hope that such a term can describe also 
the neighborhood of periodic orbits. This might include effects correspond- 
ing to longer periodic orbits close to reference orbit p. One can hope that 
a few short periodic orbits with sufficient number of correction terms can 
accurately predict energies and resonances. Numerical calculations using the 
method outlined in the paper are in progress. 
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Figure 1: The periodic orbit in the potential U{y) = — |A^y^ — ^oty^. 



Figure 2: The bilhard system of the second example. 
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